QUASI-LINEAR STOKES PHENOMENON FOR THE 
HASTINGS-MCLEOD SOLUTION OF THE SECOND PAINLEVE 

EQUATION 

A.A. KAPAEV 

Abstract. Using the Riemann-Hilbert approach, we expHcitly construct the 
asymptotic SP-function corresponding to the solution y ~ ± y^— x'/2 as ^ oo 
to the second Painleve equation y^x = 2i/^ + xy — a. We precisely describe 
the exponentially small jump in the dominant solution and the coefficient 
asymptotics in its power-like expansion. 



1. Introduction 
The second Painleve equation, 
(P2) Vxx = 2y^ + xy - a, a = const, 

was introduced more than a century ago in a classification of the second order ODEs 
Vxx — R{x, y, yx) with the Painleve property pp. To this date, equation P2 has found 
interesting and important applications in the modern theory of nonlinear waves 
[311], plasma physics 4 , bifurcation theory 5 , random matrices and combinatorics 
IHlEj, theory of semi-classical orthogonal polynomials 8 and others. 

Among its various solutions, we distinguish those with a monotonic asymptotic 
behavior as a; — > ±00 (look for complete list of asymptotics and connection formulae 
to transcendent solutions of P2 in El)- For instance, if a = and x —t +00, there 
exists a 1-parameter family of solutions approximated with exponential accuracy 
by decreasing solutions of the Airy equation, yxx = xy, see [2], 

(1.1) y^aki{x) = ^x'^'^e~"r-"\i + 0{x-^'^)), a = const. 

If in the above asymptotics < a < 1, then the Painleve function is bounded and 
oscillates as x —f —00, 

y ~ 6(-x)-i/4sin(|(-x)3/2 - ^ 

with real constants b and (j) determined by a, cf. 0. If a > 1, then the solution 
blows up at a finite point. If a = 1, then the solution of P2{oi = 0) grows like a 
square root as x ^ —00, y ~ ^J—x/2, see 

For the first time, the asymptotic as \x\ — > 00 behavior of the Painleve tran- 
scendents in the complex domain was studied by Boutroux PO]- Generically, the 
Painleve asymptotics within the sectors arga; e (|n, |(n -I- 1)), n € Z is described 
by a modulated elliptic sine which trigonometrically degenerates along the direc- 
tions argx = ^n. 

Besides generic solutions, Boutroux described 1- and 0-parameter reductions of 
the trigonometric asymptotic solutions which admit analytic continuation from the 
ray arga; = into an adjacent complex sector. As x 00 in the interior of 
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the relevant complex sector, such solution is represented in the leading order by a 
power sum of x and of a trans-series which is a sum of exponentially small terms, 

(1.2) y ~ (power series) + (exponential terms). 

Expansions of such kind can be obtained using a conventional perturbation analysis. 

In [T^ , it was observed that asymptotic solutions of the form (|1.2|l exhibit a 
quasi-linear Stokes phenomenon, i.e. a discontinuity in a minor term with respect to 
argx. The first rigorous study of the quasi-linear Stokes phenomenon associated to 
solutions y^a/x as x^oo for P2 and y ~ \J —xj^ as a; ^ 00 for P\ is presented 
in I .'i 1 1 . where the reader can find a discussion of other approaches to the same 
problem. 

Below, we study a quasi-linear Stokes phenomenon for the second Painleve tran- 
scendent with the monotonic asymptotic behavior y ^ \J —xj2 as x — *■ —00. Our 
main tool is the isomonodromy deformation method El ^| in the form of the 
Riemann-Hilbert problem approach via the nonlinear steepest descent method fT51 . 

As in [13. 114) . we pursue a two- fold goal, i.e. (a) we bring an exact meaning to 
the formal expression H1.2|) , and (b) we evaluate the asymptotics of the coefficients 
of the leading power series in 11.2|l . 

2. Riemann-Hilbert problem for Pi 

The inverse problem method in the form of the Riemann-Hilbert (RH) problem 
was first appHed to an asymptotic study of Pi in Further study of the RH 

problem can be found in [THIITOI 1^ IHl 1^ . 

According to ^B] , the set of generic Painleve functions is parameterized by two of 
the Stokes multipliers of the associated linear system denoted below by the symbols 
Sfc, A; = 0, 1, 2. As it is shown in [221, so(l + sosi) ^ then the asymptotic solution 
of the RH problem within the sector arga: e (^, tt) can be expressed in terms of 
elliptic ^-functions which in turn yields an elliptic asymptotics to the Painleve 
function itself. Assuming that the condition sq = holds true, we arrive to an RH 
problem which leads to a decreasing asymptotics y ^ a/x, see |13| . 

In the present paper, we first construct an asymptotic solution to the relevant 
RH problem as |a;| 00, argx £ (^, tt) assuming that 

(2.1) l-Hsosi=0. 

Below, we use the RH problem of Flaschka and Newell ^1 modified as it is 
proposed in T^. This RH problem comes in a standard way from a collection of 
properly normalized solutions of the Lax pair for P2, 

(2.2a) ' " ^ ^ ^ ^^''''^^ " ^^^^ +j)^2- Zy.ai, 

(2.2b) —^,-^ = -^Xa:,-ya2, 

where ^3 = ( J -"i ) , ^2 - ( ° V ) ' ^1 = ( ? J ) ■ 

Let us introduce the piece- wise oriented contour 7 = CUp+ U p_ U|^q 7^ which is 
the union of the rays 7^ = {A g C: |A| > r, arg A = f + f (^ — 1)}, fc = 0, 1, . . . , 5, 
oriented toward infinity, the clock- wise oriented circle C = {A€C: |A|=r}, and 
of two vertical radiuses p+ = {A € C: |A| < r, arg A = ^} and p- — {X £ <C: |A| < 
r, arg A = — f } oriented to the origin. The contour 7 divides the complex A-plane 
into 8 regions fife, k £ {left, right, 0, 1, 5}. r^ioft and ilnght are left and right 
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halves of the interior of the circle C deprived the radiuses p+, p-. The regions Qk, 
fc = 0, 1, . . . , 5, are the sectors between the rays and jk-i outside the circle, see 
Figure O 




Figure 2.1. A Riemann-Hilbert problem graph for 4" (A) associ- 
ated to P2- 



Let each of the regions ilk, k — right, 0, 1, 2, be a domain for a holomorphic 2x2 
matrix function ^'^(A). Denote the collection of 5'fe(A) by ^'(A), 



(2.3) 



*(A)| 



xen^ 



= *fe(A), *(ci"A) =CT2*(A)a2. 



Let 5*+ (A) and ^'-(A) be the limits of 'I'(A) on 7 to the left and to the right, 
respectively. Let us also introduce two matrices ct+ = (00)' ~ (lo)- 
The RH problem we talk about is the following one: 

(1) Find a piece-wise holomorphic 2x2 matrix function ^'(A) such that 

(2.4) ^'(A)e^''^ ^ /, A ^00, = i{j\^ +xX), 
and 

(2.5) ||*Hght(A)A-""^ II < const, A ^ 0; 

(2) on the contour 7, the jump condition holds 

(2.6) ^+{X) = ^-{X)S{X), 

where the piece-wise constant matrix S{X) is given by equations: 



on the rays 7^ , 



(2.7) 5(A) 



7fc 



Sk, 



5 



2fc-l 



S2k-icr- 



with the constants Sk satisfying the constraints 
(2.8) Sk+3 = -Sk, si - S2 + S3 + S1S2S3 = -2 siuTra; 
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- on the radiuses p±, S{X) is specified by the equations 

(2.9) Xep-: 4',eft(e'"*A) = *Hght(A)M, 

A e p+ : l'Hght(A) 1'ieft(A)cr2MCT2, 

where 

(2.10) M = (e*-("-5)-3 + j^^^^ 
with 

J+ = if i + a^N, i.e. a ^ { i, |, . . . } , 
J_=0 if i-a^N, i.e. a ^ {-i, -f , . . . } ; 

- on the circle C, the piece- wise constant jump matrix 5(A) is defined 
by the equations 

(2.11) 

*o(A) = ^'HgM(A)£;%\ *i(A) - *Hght(A)i;, *2(A) - *Hght(A)i;5i, 

*3(A) = ^'ieft(A)a2i;cr2 5'3"\ 4'4(A) *ieft(A)cr2£^cr2, 

^'5(A) = ^-irft (A) 0-2^^(72 5*4. 

The connection matrix i?, the Stokes matrices 5*^ and the monodromy 
matrix M satisfy the cychc relation, 

(2.12) ES1S2S3 = (J2M-^Ea2. 

A solution of the RH problem H2.4|l - H2.12|l . if exists, is unique. 

Remark 2.1. For the A-equation associated with PII (|2.2al) as well as for the 
above RH problem, the values a = ^ + n, n g Z, are called resonant since the 
corresponding ^'-function may have a logarithmic singularity at the origin. It is a 
quite common fallacy, that the RH problem is not uniquely solvable for resonant 
values of such parameters. As a matter of facts, monodromy data form a locally 
smooth complex surface (12. 8|) with the special points a = i + n, n G Z, si = 
= S3 = (—1)"+^. To the latter, one has to attach a copy of CP^. A complex 
parameter describing the attached space CP^ can be interpreted as the ratio of a 
column entries of the connection matrix E 16 . Thus neglecting the attached space 
CP may lead to the loss of uniqueness in the RH problem solution. 

The asymptotics of ^(A) as A ^ 00 is given by 

(2.13) y(A) *(A)e^"-' ^ I + (-iJCag + |(Ti) + ^(A-^), 
where 

(2.14) 5^ = - - W + ay 

is the Hamiltonian for the second Painlcvc equation. Thus y{x) can be extracted 
from the "residue" of Y{\) at infinity, 

(2.15) y = 2 lim AYi2(A) = 2 lim Ay2i(A). 

A — >oo A — >oo 

Equation H2.15(l specifies the Painleve transcendent as the function y = f{x, a, {sk}) 
of the independent variable x, of the parameter a in the equation and of the Stokes 
multipliers Sk (generically, the connection matrix E can be expressed via Sk using 
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equation H2.12|l modulo arbitrary left diagonal (or triangular for half-integer a) 
multiplier; at the special points of the monodromy surface (|2.8|l si — ~S2 — S3 — 
(— 1)"+^, a = ^ + n, n e Z, the connection matrix E contains a parameter r E 
CP^ which specifies a relevant classical solution of ^2)- Using the solution y = 
f{x, a, {sk}) and the symmetries of the Stokes multipliers described in | 23) . we 
obtain further solutions of P2: 

(2.16) y ^ -f{x,-a,{-Sk}), y = f{x,a,{s^}), 

y = ei^"/(ei^"x,a,{sfc.,24), 
where the bar means the complex conjugation. 

3. RiEMANN-HlLBERT PROBLEM FOR 1 + SqSi = 

First of all observe that our RH problem can be transformed to an equivalent 
RH problem with the jump contour consisting of three branches, see Figure 13.11 
(the method of transformation of Riemann-Hilbert graphs is explained in detail in 




Figure 3.1. A transformed RH problem graph. 



Let us assume now that 

(3.1) l + SQSi-O. 
Constraints (|2.8|) imply that si ~ sq = — 2sin7ra as well, so that 

(3.2) si = e-''"'^°'+i\ So = -l/si = e"'-("-5), = 1, 
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while the parameter S2 remains arbitrary. Using , it is straightforward to check 
that the jump matrices in the transformed jump graph in Figure mi are as follows, 



(3.3) iSiSoSi)-' = 

SiSoSiS20'2Si ^0-2 = 

a2SiSoSiS2fJ2S^ ^ = 





-si 


Sl Sl(si+S2 

Sl(si+S2) -Sl 


For the connection matrix E we have: 



^ I , 0'2(5iS'oS'i) ^(T2 = 



Sl 
-1/si 





Sl 



1/si 



-1/si 


-Sl 





1 Sl + S2 

1 
1 

-Sl - S2 1 



(3.4a) a - i ^ Z, 

■=-1, Sl 
■ = +1, Sl 



Si(si+S2) 1 I 
. ( 1 



where p is constant; 

(3.4b) a-^^neZ, si = -sq = (-l)"+\ Micr2 = (-1)"-?^ + J±ct± : 



1 0^ 
1 h 



n e Z>o, J+ 7^ 0, J_ = 0: ^^^o^^^j^i = p''H(j2 

where p and q are constants. 

Ifa— i=rieZ and si = — sq = (—1)"+^ then the condition J+ = J_ = is 
equivalent to the equation S2 = (—1)". On the one hand, these points have impor- 
tant geometrical meaning as being special points of the monodromy surface H2.8|l . 
On the other hand, in this case, equation (|2.12|l does not provide any restriction to 
the connection matrix E which becomes arbitrary. Taking into account that, for 
a — ^ = n G Z, the matrix E is determined up to an upper or lower in dependence 
on the sign of n triangular left multiplier, it contains one essential parameter which 
parameterizes a family of the (classical) Painleve functions. 



4. Case a - i ^ Z and 1 + sqSi = 

In the case of the non-special point of the monodromy surface, even for a half- 
integer a, the entries of the connection matrix E do not affect the Painleve function, 
and thus the jump graph can be deformed to one depicted in Figure ITTl Here, 

(4.1, ^^{x "„'■). 

and other jump matrices are as above. For our convenience, we put the node points 
of the jump graph to the points Ai,2 = iy—aT^j arg a; G [^j'"']. It is worth to note 
that the jump graph consists of the level lines \m.g{\) = const, Heg^X) = const 
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Figure 4.1. An RH problem graph for non-special points of the 
monodromy surface. 



emanating from the critical points A = —xjl and A = for the function 

(4.2) 5(A)=z|(A2 + |)^/^ 

The function ^(A) is normalized at infinity by the asymptotic condition H2.4|l . 
Since we eliminate from the jump graph the circle around the origin, the asymptotics 
for 5* (A) at the origin (|2.5I) has to be replaced by the condition 

(4.3) ||*(A)S'iS'o£;"^A""''^' II < const, A ^ +0. 

The RH problem depends on the free parameter S2 due to a jump across the 
imaginary axis. Using a quadratic change A^ = ^ supplemented by a gauge trans- 
formation of the ^'-function, it is possible to obtain a disjoint jump graph. Unfor- 
tunately, unlike the cases considered in |13[I14| . in this case, one meets a difficulty 
with the normalization of the transformed RH problem which results in ambiguity 
in the asymptotics of the Painleve transcendent. Thus we prefer to deal with the 
connected jump graph shown in Figure EtI 

4.1. Reduced RH problem with 1 + sqSi = si + S2 = 0. Consider the non- 
special RH problem for P2 corresponding to 1 + sqSi = si + S2 =0. Then the 
non-specialty assumption implies that a — ^ is not integer. The RH problem jump 
graph coincides with one depicted in Figure ini except for the jump across the curve 
lines emanating from the origin and approaching the vertical direction since now 

S^S^^ = 02828^^02 — I- 

The reduced RH problem is formulated as follows: find a piece- wise holomorphic 
function 4* (A) such that 

(1) ^(A)e^'^=' -> /, A -> 00; 
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(2) ||l'_(A)5iS'o£^o ^""'^'11 < const, A ^ 0, where So is the connection 
matrix E defined in H3.4|l corresponding to si + S2 = 0; 

(3) ^(A) is discontinuous across the contour shown in Figure UTI fwith the 
trivial jumps across the hues emanating from the origin and approaching 
the vertical directions). 

Theorem 4.1. Let aigx e [^,tt], a — ^ ^ Z and 1 + sqSi = si + S2 = 0. Then, 
for large enough \x\, there exists a unique solution to the reduced RH problem. 
If additionally the parameter a G {+1,-1} is defined in such a way that si = 
^-i-iTij{a+^) ^/jg^ ^/jg corresponding Painleve function has the asymptotics 

(4.4) y{x) = 2/1 (x, a, a) = a^e--f + 0{x~^'^), ^ oo, argx G , tt] . 

Proof. Uniqueness. Let the reduced RH problem admit two solutions ^'i(A) and 
^2('^)- Since all the jump matrices are unimodular, their determinants det^j(A) 
are continues across the jump contour, bounded at the origin and therefor are entire 
functions. Using the Liouville theorem and the normalization of ^'(A) we conclude 
that det5'j(A) = 1. Therefore there exists the ratio x(A) = >l'i(A)>I'^^(A) which 
is continuous across the jump contour for 'I'j(A), remains bounded as A and 
is normalized to the unit matrix as A — > oo. Therefore, by the Liouville theorem, 
X(A) ^ /. 

Existence. Consider an RH problem on the curve line segment [—^J—x/2, 0) U 
(0, yJ—x/2], where arg yJ—x/2 e [^f j O] : see Figure lOl (because ^t^-function can be 
analytically continued to C\{0}, the jump contours in the RH problem can be bent 
in any convenient way and even "kiss" the infinity), with the quasi-permutation 
jump matrices: 




Figure 4.2. Jump graph to the RH model problem 1. 
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RH model problem 1. 

(1) $"(A)e'''"^' ^ /, A^c5o; 

(2) ||$0(A)S'iS'o£;(7^A-"'"-^ II < const, A ^ 0; 
(3) 



$° (A) = $° (A)a2Sa2, A G (-^-2^/2,0), 
<i>0(A) = <i>"_(A)E-i, Ae(0,V^), 



oy ' \si 

A solution to this RH problem is found explicitly in [21] in some different no- 
tations. Namely, on the complex A-plane cut along [— a/— x/2, —x/2], define the 
following scalar and matrix functions I3,Y,5: 



(4.5) /3(A) 



A - V'^^\i/4 



whose branch is fixed by the condition 

/3 ^ 1 as A ^ c 



(4.6) 



era 



(4.7) SiX)^( ^^-'^ ^y, . = --Lln(...0= 



where the branch of the square root is fixed by its asymptotics y/X^ + ^ = A + 
0{\^^) as A — > +00, and z'^ is defined on the plane cut along the negative part of 
the real axis and its main branch is chosen. 
As easy to see, 

(4.8) S{~\)S{\) = 1. 

Furthermore, the functions introduced above enjoy the following jump properties: 



(4.9) y+(A) =r_(A)(-ia)ai, A e (-7=^, 7=^), 

(4.10) <5+(A)<5-(A) - e2^■'^ A G (-^^^,0), 
5+(A)<5_(A) = ^ ^ (0, v/^^), 



(4.11) 5+(A) = -9_(A), Ae(-v/^,vW2). 
Therefore the function 

(4.12) $° = YS^^'e-s"'' 
satisfies the jump conditions of model RH problem 1. 
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Because of the asymptotics as A ^ oo 

1 



(5(A) = l + 0(A-2), 

(4.13) giX)^i{p^ + x\)+i^+Oi\-^), 

we find the asymptotics as A — > oo of $"(A), 

(4.14) $"(A) = (/ + ^{-^<J3 + (T^^^i) + 0{X'' 
Using asymptotics as A 0, 

/3_(A) = e-''"/4 + 0(A), 

r-(A) = -^(/ + *aai+0(A)), 

5_(A) ^ e-''"'{~2xyX-^''{l + 0(A2)), 

(4.15) ff_(A) = -^{-xf/^ + V2{~x)^/^X^ + 0(A4.t-i/2), 

it is easy to see that 

(4.16) $°_(A) = C(/ + 0(A))A-'""'"^ 

with some constant matrix C, therefore, using iTOa . \\<i>"_{X)SiSoEo^X-'^''^\\ < 
const. Thus the function <i>°(A) H4.12|l solves model RH problem 1. 

Remark 4.1. A similar quasi-permutation RH problem in an elliptic case is solved 
in 22 . More general quasi-permutation RH problem is solved in |24| . 

For the subsequent discussion it is also worth to note that 

(4.17) CT2$"(-A)cr2 $"(A). 

We also point out that <i>°(A) is singular at A = ±y/—x/2. To "smoothen" this 
singularity, let us consider the following 

RH model problem 2. Find a piece-wise holomorphic function $(A) with the 
jumps indicated in Figure IT!^ Since we do not normalize *i>(A), it is determined up 
to a left multiplication in an entire matrix function. If we will use this solution in 
a domain different from C, then <E'(A) is determined up to a left multiplication in a 
matrix function holomorphic in this domain. 

It is possible to construct a solution ^(A) of this model problem using the classical 
Airy functions. The Airy function Ai(z) can be defined using the Taylor expansion 

EH EH], 

(418) Ai(.)-^— y^^£i^i±il^ l^yS^Tik+l)^ 

^ ' ^' 3V3r(|)^4^^ r(i)(3fc)! 3i/3r(i)^4^^ r(|)(3fc-fi)! • 
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Figure 4.3. RH model problem near A = -x/2. 

Asymptotics at infinity of this function and its derivative are as follows, 
(4.19) Ai(.) . ^e~l^'''{Y^yir,-^-^-^^z-^ 

n— \2 f 

as 2 — > oo, arg z e (— vr, 7r). 

Introduce the matrix function Zo(A), 
where 

(4.21) V2{z) = e*2"/3Ai(e^2V3^)^ ^^(^) ^ ^,(^)^ 
As |z| oo, argz G tt, ^), this function has the asymptotics 

(4.22) Zo(z) =z--/4i= (J _\) (/ + 0(z-3/2))ei^^^^-. 
Also, introduce the matrix functions 

(4.23) Zi(z) Zo(z)Go, Z^^z) Zi(z)Gi, ^3(2) ^2(^)62, 

Go = G2 - (^_^ 5) , Gi = (^J . 

By the properties of the Airy functions 123123, Zj[z^ have the asymptotics (I4.22() 
as z — > 00 and arg z e (— tt + ^j, ^ + ^j) , j = 0, 1, 2, 3. 
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Define a piece- wise holomorphic function Z (z), 

(Zoiz){tsir^/^ arg^e(-f,0), 

(4.24) Z^"{z)= lz,{z){is^r^/^ argze (-f + f j,f j), J = 1,2, 

[Z,{z){zs,r^/^, argze(f,f). 

By construction, Z^^{z) has the uniform asymptotics 

(4.25) Z«^(z) = ^""'/'^ (\ (/ + 0(z-3/2)^g|.^/V3 (,5^)^3/2^ ^ _ 
and jump properties indicated in Figure ^31 



oo, 




Figure 4.4. A model RH problem solvable by the Airy functions. 



Using the change of the independent variable 



(4.26) 

and the notation 
(4.27) 



we compute the "ratio" of $° and Z^^ in the annulus ci\x\ 2+'^ < |^| < C2|a;| 2+*^^ 
where < ci < C2 and j < ^ < ^ a-re some constants, 

(4.28) $0(A)(Z^^(z))"' - (/ + ©(x-i/^C^/^) + 0(^-1/2^-2)) ^ 

X (I \) (i^.., + i±^/)e^f-2i-(-./2)-/^ 

Define the matrix function $,.(A), 
(4.29) 

<i>„(A) = ^ (I + i±^/)e^f-2§-(-x/2)-/4z«^(z(A)). 
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Let US introduce the piece- wise holoniorphic function $(A), 







r$,(A), 




-x/2 < R, 


(4.30) 


*(A) = 1 


1 a2'I>r(e" 


^"A)a2, 


\X+y^-x/2\ < R, 






[$°(A), 


|A± V- 


-x/2\ > R, 



where R — c\x\ 2"'"'^ for a constant c > and j < ^ < \- We look for the exact 
sohition of the reduced RH problem in the form of the product 

(4.31) *(A) = x(A)<i>(A). 

The correction function x(A) solves the RH problem whose jump graph is shown in 
Figure ^31 The jump matrix across the circle centered at A = is given by 




Figure 4.5. RH problem jump graph for the correction function x(A)- 



the ratio $°(A)<i>~^(A) which, due to (|4.28|) . I4.29|l and choosing e = 2/5, satisfies 
the estimate 



(4.32) ||$"(A)$7i(A)-/|| <c|x| 



-3/10 



The jump matrices across the exterior parts of infinite contours approach the unit 
matrices exponentially fast, 

(4.33) ||$°(A)5,-($°(A))"' - /|| < C|x|3/20e-^l-l'''l«l'''. 

Now, the solvability of the reduced RH problem is straightforward. Indeed, 
consider the system of singular integral equations X- = ^ + Kx~ equivalent to the 
RH problem for x(A), 

(4.34) ^_^X)=I+^ j iO {GiO '^^ 



2TTi J^^ ^^'^ ^S-A- 
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where 7 is the contour shown in Figure ESI and G(A) = $_(A)5j$+^(A) is the 
relevant jump matrix. Using the estimates (|4.32|l . (|4.33() and the boundedness 
of the Cauchy operator in L^(7), the singular integral operator K, which is a 
superposition of the operator of the right multiplication in G{\) — I and of the 
Cauchy operator C_, satisfies the estimate, ||^||^2(..y) ^ c\x\^'^^^ , where the precise 
value of the positive constant c is not important for us. Thus K is a contracting 
operator in (7) for large enough |a;| . Since KI is a square integrable function on 7, 
and observing that p — X- ~I satisfies the singular integral equation p = KI + Kp, 
we find p G ^^(7) and therefore X- = I + P iterations. 

The solution found by iterations implies the asymptotics of x(A) as A — > 00, 

(4.35) x(A) = /+y(^(x-2/5) + 0(A-2). 

A 

Using H4.35|l and (|4.14|) in H4.31|) and taking into account the expansion H2.13|l . we 
see that the Painleve function corresponding to the reduced RH problem has the 
asymptotics as |a;| 00, argx € [^,7r], 



(4.36) y{x) = e-"| + 0{x'^/^). 

This completes the proof. □ 

4.2. RH problem with 1 + SqSi = and Si + S2 7^ 0. Let us go to the case of 
1 + sqSi = and arbitrary si + S2 ^ 0, see Figure HjI We look for the solution 
^'(A) in the form of the product 

(4.37) *(A)=x(A)*(A), 

where ^(A) is the solution of the reduced RH problem, i.e. with Si + S2 = 0. 

The correction function x(A) satisfies the RH problem on two level lines Im g{X) = 
const emanating from the origin and approaching the vertical direction, see Fig- 
ure ^21 

(1) x(A)^/, A ^00; 

(2) ||x(A)*(A)S'i5'o-B-iA-"'"3 1| < const, A ^ +0; 

(3) x+(A)=X-(A)G(A), 

G(A) ■.= Gi{X) = ^-{X)S2S^^^z\X), Xe (0,+ioo), 
G(A) := G2(A) = ^_(A)a25254-V2#li(A), A e (0, -too), 
a2Gi(-A)(T2 = G2(A). 

Theorem 4.2. If a — i ^ Z, 1 + sqSi = 0, argx e [^,7r] and \x\ is large enough, 
then the RH problem jjj)) is uniquely solvable. The asymptotics of the 

corresponding Painleve function as |x| —> 00 m the indicated sector is given by 

(4.38) y ^ yi{x,a,cr)- 



where yi{x, a, a) ^ e^^"^ x/2 is the solution of the Painleve equation correspond- 
ing to \ + spsi = si + S2 = 0, while the parameter a G {+1, —1} is defined by the 
use of the equation si ~ e^^^'^i'^+s) _ 
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Figure 4.6. An RH problem graph for the correction function x(A). 

Proof. It is enough to prove the solvability of the RH problem for x{^) above. In 
a neighborhood of the imaginary axis, ^(A) — x(A)<I>"(A), see (|4.31|) and H4.30|l . 
Hence, using H4.12|l . we have the following expression for the jump matrix across 
(0, +ioo), 

(4.39a) Gi(A) = / + (si + S2)5'^e-^3^YG+Y-\-^ , 

(4.39b) G2(A) = / - (si + S2)5-^e^3xY^-Y'^x~^- 

Because the jump contour coincides with the level line lu\g{\) ~ const, the jump 
matrix Gi(A) exponentially approaches the unit matrix as A — > +ioo. Similarly, 
G2(A) — / decreases exponentially as A — > —ioo. For the exponential e^*^'^-', due to 
(|4.15|l . the origin is the saddle point. 

The jump matrices have algebraic singularity at A = 0, 

Gi(A)=/+(si+S2)e-'"^"(-2a;)-2'^A4'^e-'^(-^)'''+2^(-")'''^'x 
(4.40a) X i(zaa3 + fTi)(l + 0(A) +0(2;-9/i")), 

G2(A) - (si +S2)e'^"''(-2a;)-2^A4^e-¥(--)'''+2V^(-)'''^')x 
(4.40b) X i(iaa3 + CTi)(l + 0(A) + 0(a;^^/i")). 

Consider the model RH problem with the jump matrices Gj(A), j = 1, 2, of the 
form H4.39|l but with the matrix x(A)y(A) replaced by its asymptotics at A = 0, 
i.e. 

(4.41a) Gi(A) = / + (si + S2)S^e-^9 + ^^)^ 

(4.41b) G2(A) = / - (si + S2)'5-'e2f i(iaa3 + di). 
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Assume for a moment that 
(4.42) Rev>0, i.e. Re(CTa)>0. 

Since Gj{X) — / are constant nilpotent matrices multiplied in scalar functions, so- 
lution of the jump problem is given by the Cauchy integral, 

(4.43) 

The condition (|4.42|l ensures the convergence of the above integrals. Obviously, 
x(A) ^ / as A — > oo. As to asymptotics at the origin, the same condition (|4.42|) 
ensures that ||x(A)|| < const, and \\x{X)'^{\)SiSoE-^ X-"""^]] < const. 
We will look for x{X) in the form of the product 

(4.44) x(A) - X{X)x{X). 

For the correction function ^(A), we have the RH problem with the jump contour 
shown in Figure but with different jump matrices: 

(1) X{X)^I, A^oo; 

(2) ||^(A)|| < const, A +0; 

(3) X+iX) - X_(A)i/(A), HiX) = x-(A)G(A)G-i(A)x='(A). 
In the equivalent singular integral equation, 

(«5) .Y.(A).7 + J-|^x_«)(//(0-/)j^, 

where 7 is the jump graph in Figure 14. 61 and A_ is the right limit of A on the 
jump contour, or, in the symbolic form, X_ = / + KX^, the operator K is the 
superposition of the right multiplication in iJ — / and of the Cauchy operator C_ . 
Using the boundedness of C_ in ^^(7), we estimate the norm 
(4.46) 

< c\\H - /||^,(^^ < c'|x|-(5-+l)/2e-f l-l^''™-^ i(argx-.) < c"\x\''^^, 

where c, c' and c" are positive constants whose precise value is not important for us. 
Taking into account the assumption Re > (|4.42|) , the singular integral operator 
K is contracting for large enough |a:;|, argx £ tt] , and therefore equation (|4.45() 
is solvable by iterations. 

To find the asymptotics of the relevant Painlevc function, it is enough to find 
asymptotics of x(A) H4.43|l as A ^ 00, 

(4.47) x(A) = / - '^-^ 5^e-^^ dC{l + 0{X-^)) [laa^ + a,) = 

2vriA Jo 

= / + X~^h{iaa-i + cTi) + 0(A-3), 

h = -fi±fa2-5-5r(2^ + i)(-.x)-3-le-^(--)^'^(i + 0(x-i/2)). 
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Thus the asymptotics of 5" (A) as A ^ oo H2.13|l is given by 

(4.48) *(A)e''''^ = Xx'^e'^"' =1 + \-^{-i'Ka3 + |cti) + ©(A^^) = 

where hj — ft,(l-|-0(x~^''^^/"')), j = 3, 1, involves the contribution of the correction 
function X{\). Comparison of two hnes in (|4.48|) yields the asymptotics of the 
Painleve function, y = yi+ 2hi, which turns into H4.38|l for Re (era) > 0, a — ^ ^ Z 
substituting v = aa/2. 

Validity of the asymptotic formula (|4.38|l without the restriction (|4.42|) follows 
from the observation that H4.38|l is invariant with respect to Backlund transforma- 
tions. Indeed, the Schlesinger transformation of the ^'-function, 

(4.49) *(A) = i?(A)'J'(A), 

2jA yx - e{y^ + ^ ) 

yields the \l/-function, associated to a new Painleve transcendent 

(4.50) y = y + eq^, ijx ^ Vx + <^{y^ - v'^), a = -a - e, 

characterized however by the same Stokes multipliers. The latter with the def- 
inition of the parameter a via si = e^*''°"("+2) = si = e^"^'^("+2) implies the 
transformation of the parameter a, a = —a for either value of e = ±1. If both 
Re (era) > and Re {aa) = Re (era) + ere > 0, then the asymptotics of both y and 
y is described by (|4.38|l (supplemented by tilde over y, a and a if necessary). This 
observation constitutes the invariance of the asymptotic formula with respect to the 
Backlund transformation. Since the latter transformation is a bi-rational transfor- 
mation in the space with coordinates {y, yx), the above invariance can be confirmed 
algebraically. Therefore, if Re {aa) < while Re (era) > 0, the asymptotics of y 
and y are described by the formula (|4.38|) (supplemented by tilde if necessary) . The 
iterated use of the Backlund transformations completes the proof. □ 

4.3. The increasing degenerate Painleve functions. Applying the second of 
the symmetries (12.1611 to (|4.38|l and changing the argument of x in 2tt, we obtain 

Theorem 4.3. If a — ^ ^ Z, 1 + sqSi = 0, argx £ [tt, ^] and \x\ is large enough, 
then the asymptotics of the Painleve function is given by 

(4.51) y^yo{x,a,a)- 

_ £i_iO .a-|p( 1 ^ 1 ^_2^(e--.)3/^ ^ 0(X- 1/4)) ^ 

where yQ{x,a,a) = yi{x,a,a) ^ e^™ x/2 is the solution of the Painleve equa- 
tion corresponding to 1 + sqSi = S2 — sq — 0, while the parameter er G { + 1, —1} is 
defined by the use of the equation Si ~ £^"'^("+2). 

The solutions yi{x, a, cr) and yo{x, a, cr) = yi{x, a, er) are meromorphic functions 
of a; € C and thus can be continued beyond the sectors indicated in Theorems 14.21 
and 14.31 To find the asymptotics of the solution yQ[x,a,a) in the interior of the 
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sector arga; G [^,7r], we apply (|4.38|) with 52 = 50. Similarly, we find the asymp- 
totics of the solution yi{x, a, a) in the interior of the sector argx e [tt, ^] using 
(|4.51|l with S2 = — si. Either expression implies that, if oo, a.rgx e 

(4.52) yQ{x,a,(T) - yi{x,a,(T) ^ 

where we have used the relation 

, s 2i7rcr 

(so + si)r(2 + aa) = --jj ^. 

f (2 - craj 

The formula (|4.52|l constitutes the quasi-linear Stokes phenomenon for the increas- 
ing degenerate asymptotic solution of P2. 

Due to exponential decay of the difference (|4.52() in the interior of the sector 
argx e ^), we have solutions of P2, 



(4.53) y = yi(a;,a,a) ~cr\/e-"2:/2, ^ 00, argxe[^,^), 

for So = -e*'''"(°+^', si g-^'^'^C^+D, S2 = -e-"^'"("+^\ 

and 



(4.54) y = yo(a;,a,cr) ~crv/e"*'^2;/2, ^ 00, arg2;e(^,^], 

for So = ~e''"''^°'+i\ si = e-''"^^°'+^\ S2 = -e''''^("+^). 
Applying the rotational symmetry of (|2.16() to yi and yo, we find solutions 



(4.55) y = y2n+i{x, a, a) := e'^"yi(e^^"a;, a, a) ~ a(-l)" v/e-"a;/2, 
|a;|^oo, argx e - ^n, ^ - ^n), 
for S2„ = -e'^"("+^\ S2„+i =e""'^("+^), S2„+2 = -e''""("+^), 

and 



(4.56) y = y2n{x, a, a) e'^"yo(e'^"a;, a, a) ~ a(-l)'Ve~"a;/2, 

\x\^oo, argx e - ^n, ^ - ^n], 
for S2„ = -e^""("+^\ S2„+i = e-'"'^("+^), S2„+2 = -e'""("+^). 
Comparing the Stokes multipliers, we observe the symmetries 

(4.57) j/„(x, a, a) = y„+6ix, a, cr), y2n-i{x, a, a) = y2n{x, a, -a). 
Equation (|4.52l) and definitions 14.55|l and (|4.56() imply the differences take place, 

(4.58) i/2n(a;, a, cr) - y2n+i(a;, a, cr) = 



T{\ - aa) 

which are exponentially small in the interior of the indicated sectors. 



x| arga; e - ^n, ^ - ^n] 
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Using the second of the equations l|4.57|l . we estabhsh the analytic continuation 
of the asymptotics as \x\ — > oo, 



(4.59) y2n-iix,a, = j/2„(x, a, (-1)"(t) ~ a^e---x/2, 

\x\~^oo, argxe (^(l-n),^(3-n)). 

5. Coefficient asymptotics for a — | ^ Z 

An elementary investigation shows that the equation P27 Uxx = 2j/^ + xy — a, 
can be satisfied by the formal series depending on a and a parameter u = ±1, 

00 

(5.1) yf{x, a, a) = a^/^Y. bni^^r^"^^ + Oix-°°). 

n=0 

Given a, the series (15.111 is determined uniquely since its coefficients bn are deter- 
mined by the recurrence relation, 

(7(1 

(5.2) 60 = 1, ^1 = ^. 



n+l n+ln+2-l 



bn+2 — ^ bn — ^ bmbn+2-m — ^ ^ ^ hbmbn 

771—1 ^^1 m— 1 

Several initial terms of the expansion are given by 

„. , , , — r_, aa 1 + 6a^ a-Q!(ll -f 16a^ 

(5.3) Uf (x, a, (7 = (TV ^a:/2-^ 1 H ^ ■ ; -^r -) r= — 

73 + 708a2 4200^^ crQ;(1021 + 25040^ + 768a'') 



128(-a;)6 64^/2(-.T)l5/2 
10657 + 129918^2 + 132060a4 + 24024a'^ 
1024(-j;)9 ^ 
cra(248831 + 786304^2 + 4164000* + 49152^^) 



O 



512V2(-a;)2i/2 

To find the asymptotics of the coefficients 6„ in (|5.1|) as n ^ 00, let us construct 
a sectorial analytic function y(t), 

(5.4) y(i) - 2;2n-l(e"^^ a, (-1)""^) = 2;2n(e^"^^ a, (-1)"^), 



3 ' 3 3 



The function y{x) has a finite number of simple poles and, by construction, has the 
uniform asymptotics near infinity, y{t) ~ at/\/2 as \t\ 00. Using H5.1(l . it has the 
following formal series expansion. 



(5.5) y(t)=^J=^5„t-3". 

^ n— 

Let (a;) be a partial sum 

w-i 

(5.6) = 

^ n— 
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and f (t) be a product 

(5.7) v(^Ht) = t'^'-'m) ~ W) = ^ £ bn+Nt-'-. 

Because t'^^^^y^^'' (t) is polynomial and y{t)/t is bounded as \t\ > p, the integral 
of u^^^(i) along the circle of the radius \t\ = p containing all the pole singularities 
of y{t) satisfies the estimate 



\t\=P 



(5.8) i v'^^\t) dt < p^^-^ 6 \y{t)\dl<27Tp^^ma^\y{t)/t\^Cp 



3N 



\t\=p l*l=P 



On the other hand, inflating the sectorial arcs of the circle \t\ — p, we find that 

2 

(5.9) <f v^^^{t)dx=<f «W(t)di+ V / {v^^\t) - ti^\t)) dt. 

J\t\=p J\t\=R „=-3"'e'f "(p,fl) 

Since w(^)(t) = -^In + 0{t-^), the first of the integrals in the r.h.s. of is 

(5.10) / v^^\t)dt = TriaV2bN + 0{R-^). 

J\x\=B. 

Last six integrals in (|5.9I) are computed using (|5.4|) . H4.55|l . (|4.56() . (|4.57|l and H4.52|l . 

(5.11) ^ (vi'^Ht) - v^'^Ht)) dt = 

Jp 

+ 3(-l)^-i / ^3^^-2(yo(e"^^a,-a)-yl(e^-^^a,-a))d^ = 
r(2 - era) 7p 



r(2+cra)jp 

= za2-i^3^-^ (6-"° ~ ~ (l + OiN~'/')) + 
I r(^ - era) 

+ (^1)^6- '^^^+;" 7^ (1 + o{N-^^'))}+ 

Thus, letting R — oo, we find the asymptotics as — s- oo of the coefficient in 
(Oil . 

+ (_i)A^6-"£(^l±^eLiil(i + o(iv-i/6)) I + o(p3A^). 

r(i+cra) ^ 
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In particular, the asymptotics H5.12|l is consistent with the elementary observation 
that, for a = 0, all odd coefRcients vanish, 62^-1 = 0; for even coefficients, the 
leading order asymptotics reduces to the one-term expression, i.e. 



Remark 5.1. The error bound in H4.38|l . (|4.51|) and 14.5211 can be improved to 
Q(^x^3/2^ which implies the error bound for l|5.12|l 0{N^^). 

Remark 5.2. It is possible to prove that the asymptotic formula (|5.12ll remains 
valid for half-integer values of a as well. 
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